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ABSTRACT 

We extend our modeling of the ionization structure of magnetohydrodynamic 
(MHD) accretion-disk winds, previously applied to Seyfert galaxies, to a popu- 
lation of quasi-stellar-objects (QSOs) of much lower X-ray-to-UV flux ratios, i.e. 
smaller a ox index, motivated by UV/X-ray ionized absorbers with extremely high 
outflow velocities in UV-luminous QSOs. We demonstrate that magnetically- 
driven winds ionized by a spectrum with a ox ~ —2 can produce the charge states 
responsible for Civ and Fexxv/Fexxvi absorption in wind regions with corre- 
sponding maximum velocities of f(Clv) <^ 0.1c and u(Fexxv) <^ 0.6c (where c 
is the speed of light) and column densities Nh ~ 10 23 — 10 24 cm -2 , in general 
agreement with observations. In contrast to the conventional radiation-driven 
wind models, high-velocity flows are always present in our MHD-driven winds 
but manifest in the absorption spectra only for a ox <s —2, as larger a ox values 
ionize the wind completely out to radii too large to demonstrate the presence 
of these high velocities. We thus predict increasing velocities of these ionized 
absorbers with decreasing (steeper) ct ox , a quantity that emerges as the defining 
parameter in the kinematics of the AGN UV/X-ray absorbers. 
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merical — quasars: absorption lines — X-rays: galaxies 
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Introduction 



The launch of Chandra and XMM-Newton ushered a new era in X-ray astronomy of 
AGN outflows with the discovery of absorption lines in the spectra that enabled for the first 
time accurate charge state and velocity measurements. The long observations of a number 
of AGNs revealed transitions of charge states as diverse as Fe I through Fe xxvi. Since any 
atomic gas with bound electrons absorbs X-rays these ions span a range of ~ 10 5 in ionization 
parameteiQ £, a fact that underscores the great utility of X-ray spectroscopy. 



In a subsequent development, iHolczer. Behar fe Kaspil ( 120071 ) and iBeharl ( 120091 . here- 
after, B09) developed a statistical measure of the plethora of the transitions in the Chan- 
dra/XMM spectra, the absorption measure distribution (AMD), namely the differential 
hydrogen-equivalent column Nh of specific ions per decade of £, i.e. AMD = dNn/dlog£. 
Moreover, the AMD was found to be roughly constant, i.e. Ah to be roughly independent of 
£, in the small number of Seyferts for which the data quality allowed a quantitative analysis. 
The functional form of the AMD is significant as it can provide the plasma density along 
the observer's line of sight (LoS), which for constant AMD is n(r) oc r -1 . 



Motivated by the AMD syste matics, Fuku mura et al.l ( 20 



the photoionization code XSTAR (Kallm an fe Bautistal 12001 ) to determine the ionization 



0|, hereafter, FKCB) employed 



structure of the 2D winds of IContopoulos fe Lovelace! (11994 hereafter, CL94) which pro- 
vide for density profiles such as n oc r _1 . This density dependence on r yields also £ oc r -1 , 
thereby allowing for ionic species of decreasing ionization with distance, but of columns sim- 
ilar to those of high ionization. Importantly, these models are scale free: with the radial 
coordinate r normalized to the Schwarzschild radius, r s , and the mass flux to the Eddington 
rate, £ is independent of the black hole mass M, implying broad applicability in galactic 
and extragalactic settings. Assuming an ionizing spectral energy density (SED) of F u oc v^ 1 
between 1 and 1000 Ryd, these models were successful in reproducing the observed: (i) Slow 
velocities (v ~ 100 — 300 km s _1 ) for the low ionization transitions like Fexvil and fast 
outflows (v ~ 1,000 — 3,000 km s -1 ) for the high ionization ones such as Fexxv, and (ii) 
AMD almost independent of £ for —1 ^ log£ ^ 4, in agreement with the results of B09. 

While X-ray absorption lines in Seyfert spectra a re rather recent discoyeries, UV absorp- 



tion l ines in Seyferts and QSOs have been known (e.g 
2000). Also known since the earlier ROSAT surveys (e.g 
19961 ) is that the X-ray-to-UV flux ratio of the broad absorption line ( BAL) QSOs [i.e. QSOs 



Crenshaw et alj2003l ; lBrandt. Laor. fe Wills 
Kopko et al.lll994J : Ereen fe Mathur 



with blue absorption C IV a nd Lya troughs of Av/c ~ 0.04 — 0.1 (e.g. iHewett fe Foltzl 12003 



Srianand fe Petitjeanll2000l )] is smaller than that of the QSO majority, possibly due to absorp- 



x £ = L/(nr 2 ) where L is an ionizing luminosity (between 1 and 1000 Ryd), n is the plasma number 
density and r is distance from the ionizing source. 
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tion of the X-rays by the BAL plasma. Indeed, this was confirmed by the ASCA detection of 
high X -ray absorption column Njj > 5 x 10 23 cm -2 (iGallagher et al.lll999l ). iGallagher et al. 



(120061 . hereafter, G06) later conducted a Chandra survey combined with known UV absorp- 
tion properties that supported the earlier claims. The Chandra data of BAL QSOs indicate 
that a ox (BAL) ~ —2.21 (G06) is smaller than the mean QSO value a ox (mean) ~ —2.0 H. 
This result is augmented by a correlation between a ox and the 1 — 5 keV X-ray photon index 
T: increased photoelectric absorption of soft X-rays (i.e. smaller a ox ) also yields a smaller 
effective T, as observed. 

The high outflow velocities of the prominent UV resonance lines (Civ and Lyq ) in 
BAL QSOs were traditionally ascr ibed to radiation-driven winds ( Weymann et al. 199ll ). in 
anal ogy with the winds o f O stars (ICastor. Abbott fc Klein II1975T ) and were modeled as such 
(e.g. iMurray et al.l Il995l . hereafter, MCGV). MCGV recognized and included heuristically 
the effects of wind ionization and its shielding from the QSO X-rays, a crucial pro c ess as 



ionization reduces severely the effectiveness of line driving. iProga. Stone, fc Kallmanl (12000 



hereafter, PSK) presented 2D hydrodynamic simulations of these winds, including X-ray 
ionization, showing that the required shielding is provided by the section of the wind closest 
to the X-ray source that "failed" to launch by being too highly ionized, thereby allowing 
exterior segments to achieve velocities in agreement with C IV observations. 

However, recent X-ray observations of BAL QSOs revealed absorption features in their 
spectra identified with highly ionized Fexxv/Fexxvi of column density Nh ~ 10 23 — 10 24 
cm" 2 , blueshifted to high velocities v/c ~ 0.4 - 0.7 (e.g. APM 08279+5255, PG 1115+080 
and H 1413+117) indicating that X-ray ionization does not necessarily in hibit outflows, which 



can occur a t velocities even highe r than those seen in the UV lines (e.g., IChartas et al.ll2002 



20031 . 120071 ; IChartas et al.l 120091 . hereafter, C09). Additional X-ray studies have revealed 



a number of non-BA L QSOs that also e x hibit similar X-ray absorbers at high velocities 



v/c ~ 0.1 - 0.5 (e.g., IPounds et alJ 12003 ; iReeves et al.l 120031 : iReeves et al.l 12009 ). while in 



APM 08279+5255 C09 have also noted a correlation between T and the velocity of Fexxv. 

Motivated by these observations, we examine in this letter the conditions under which 
the magnetically-driven winds discussed in FKCB can reproduce the observed velocities of 
the BAL QSO X-ray features (Fexxv) along with those of their more common UV transitions 
(Civ). In §2 we summarize the physics of MHD accretion disk winds and the differences of 
the ionization properties between Seyferts and BAL QSOs. In §3 we present our results and 
demonstrate a number of well-defined correlations among their kinematics, column, spectral 
index, and LoS angle and we conclude with a summary and discussion in §4. 



2 The spectral index a ox = 0.384 log(/2kcv//250o) measures the X-ray-to- UV relative brightness where 

o J II 

/2koV an d /2500 are respectively 2 keV and 2500 A flux densities (jTananbaum et al.1 1979) 
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2. The MHD Disk- Wind Model 

In this section we present a brief outline of the MHD winds, originally formulated by 
Blandford & Payne (1982) and generalized by CL94 to include arbitrary distribution of axial 
current with radius. Here and in FKCB we focus on the current distribution that produces 
a density profile n(r) oc r" 1 , crucial for obtaining the observed AMD behavior. The same 
distribution leads also to a toroidal field oc r _1 that has equal magnetic energy per decade 
of (cylindrical) radius. 

Self-similarity is assumed, i.e. power-law radial dependence for all variables and solution 
of the remaining angular part of these equations. As discussed in FKCB, this assumption 
is not very restrictive and justified a posteriori by the large number of decades of £ in the 
AMD form. 

The fundamental quantity of axisymmetric MHD is the magnetic stream function \l/(r, 9), 
assumed to have the form \&(r, 9) = (r /r ) g ^/(9)^/ , with \l/ the poloidal magnetic flux 
through the fiducial wind launch radius at r = r a . *&(0) is its angular dependence to be 
solved for and q ~ 1 a free parameter that determines the radial dependence of the poloidal 
current. The scalings of the poloidal magnetic stream function carry over to the rest of the 
wind properties of which we show only the magnetic field, velocity and density (see FKCB): 

B(r,9) = (r/r o y- 2 B(0)B o , (1) 
v(r,9) = {r/r )- l ' 2 ~v{9)v , (2) 
n(r,9) ee {r/r ) 2 «-'n{9)B 2 v- 2 m- p \ (3) 

where m p is the proton mass. The dimensionless angular functions denoted by tilde must be 
obtained from the conservation equations and the solution of the Grad-Shafranov equation 
with initial values on the disk (denoted by the subscript "o") at (r = r o ,0 = 90°). The 
density normalization at (r o ,90°), setting n(90°) = 1, is given in terms of dimensionless 
mass-accretion rate rh (see FKCB) by 

2a T r s 

where r]w is the ratio of the mass-outflow rate in the wind to rh, assumed here to be of 
order unity and is the Thomson cross-section. It is important to note that because the 
mass flux in these winds depends in general on the radius, rh always refers to the mass 
flux at the innermost flow radius at r ~ r s where r s is the Schwarzschild radius. In the 
present treatment we adopt the value q = 0.93 resulting in n oc r~ L14 , the steepest density 
dependence on r implied by the AGN AMD data of B09, and in order to allow for the 
somewhat higher observed X-ray column than UV column. 

With the dimensionless, mass-invariant wind structure (see Fig. la) for given rh and 
9, the only significant difference in the wind ionization properties across objects of different 
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-6.5 log(n[cm- 3 ]) n. 5 




Equator [log(r/r sin9)] Photon Energy [eV] 

Fig. 1. — (a) Poloidal density structure log(n[cm~ 3 ]) of MHD wind with fn = 0.5; density 
contour curves (dotted lines with numbers) and magnetic field lines (solid curves) for M = 
10 9 M & . Also shown are the positions of the Civ (square = □) and Fexxv (circle = o) 
which shift outward along LoS of decreasing 9 [from 80° (innermost) to 30° (outermost) by 
a 10° increment]. Note that the Civ position for 9 = 30° lies outside the figure range, (b) 
The form of the assumed input SED consisting of a thermal MCD of innermost temperature 
kT m = 5 eV and a PL continuum of photon index T normalized by a ox . 



luminosity is the spectral distributi on of ionizing radiatio n. While in FKCB we used a 
spectrum of the form F v oc v~ l (e.g. ISim et al.ll2008l . 120101 ) . more appropriate for Seyferts, 
here we add a bright UV disk source. The spectrum used in the present work is shown in 
Figure lb; it comprises a multicolor-disk (MCD) with an innermost te mperature of 5 eV 



and a n X-ray power-law (PL) of photon index V normalized by a ox (e.g. lEverettl 120051 ; [Sim 



20051 1 . We do not include a soft X-ray excess in the SEP, a f eature more appropriate for 



narrow-line Seyfert spectra (e.g. iPounds et al.ll2003l ; ISiml 120051 ). The PL has a low energy 
cut-off at 5 eV and a high energy one at 200 keV. The total (X-ray plus UV) luminosity is 
L = 3 x 10 45 erg s -1 . 



3. Results 



With the background flow (fn = 0.5) and the spectrum of the ionizing radiation (r = 
— a ox = 2) given, we follow the same procedure as in FKCB: we split the wind logarith- 
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mically into a number of radial zones; we employ XSTAR to compute the ionization and 
opacities/emissivity in each zone along an observer's LoS. 

In Figure [2] we show the resulting distribution of the hydrogen-equivalent column den- 
sities ANh of iron and carbon for 9 = 50° as a function of £ (optimized here to model the 
outflows in APM 08279+5255) along with the corresponding LoS velocity (dashed curve), to 
be read on the right vertical axis. Note the well-defined velocity gradient of the wind with 
ionization parameter £ due to its continuous structure. The Fexxv ions identified through 
their resonance transitions of ls-2p/3p obtain their peak hydrogen-equivalent columns of 
N H (Fexxv) ~ 2 x 10 23 cm -2 at log£ > 5 with velocities v > 0.5c. On the other hand, 
the Civ ions (2s-2p transitions) yield iV#(Clv) ~ 10 23 cm -2 with corresponding LoS veloc- 
ities v ~ 0.1c, roughly consistent with the UV/X-ray absorbers in APM 08279+0255 date@ 
(C09). The obtained ionization structure directly reflects the spatial positions of these ions; 
i.e. r(Fexxv)/r s ~ 5 - 40 and r(Cw)/r a ~ 200 - 700 for 9 = 50° (see Fig.la) assuming a 
sine; le Lofl Note that the equivalent width (EW) of the modeled C IV absorption seems also 
to be fairly large due to a wide spread of its peak column distribution (i.e. 1.5 < log£ < 4) 
equivalent to Av ~ 18, 000 km s" 1 in agreement with the typical width of C IV BALs (i.e. 
Av ~ 10, 000 - 30, 000 km s" 1 ). 

1e+24 1 
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q 

< 
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Ionization Parameter log(^ [erg cm s" 1 ]) 




Fig. 2. — Simulated distribution of local column densities ANh (left ordinate) and the 
outflow velocity v (dashed curves; right ordinate) as a function of ionization parameter £ for 
carbon and iron along the 9 = 50° LoS. The shaded regions denote the parameter space (in 
£ and v) for which the local column is dominated primarily by Civ or Fexxv. 



3 It is conceivable that scatterred/reflected UV photons could externally fill in the "true" Civ absorption 
feature to seemingly reduce its intrinsic column (S. Kraemer, private communication). 

4 Diflcrent LoS for the UV/Optical and X-ray emitting regions fe.g. bai etailboid) could alternatively 
be modeled with radiation transfer in a more complex source geometry. 
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Fig. 3. — (a) Expected correlation of MHD wind velocity v with V {solid curves; left ordinate) 
and a ox {dashed curves; right ordinate) indicated by (T, a ox ) for C IV and Fe xxv with fn = 0.5 
and 9 = 50°. (b) Expected correlation of wind velocity v with ANh and 9 (indicated by 
numbers in degree) for Civ and Fexxv with V = —a ox = 2. Shared regions indicate the 
observed velocity dispersions of Civ and Fexxv for APM 08279+5255. 

The difference in the Fexxv and Civ velocities above from those in FKCB [v (Fexxv) 
~ 3,000 km s _1 and u(Crv) ~ 300 km s _1 ] begs an explanation, considering the scale 
invariance of these winds and the similar values of 0(~ 50°) presented in both cases. This 
can be traced to the different fn and/or ionizing spectrum F v . Indeed fn here is ~ 4 times 
higher than that in FKCB, leading to a large increase in jV#; furthermore, a ox is smaller here 
than in FKCB, with the p resent value a ox = — 2 consistent with the values of UV- luminous 



QSOs flSteffen et all 120051 ) 



These differences affect the Fe xxv and C IV velocities as follows: The combined increase 
in fn and decrease in a ox imply lower ionization of the plasma at the smallest r, despite the 
fact that the ratio L/n Q does not change; this is because the luminosity L in the definition of 
£ involves an integral over the spectrum, while the ionization of the gas is affected mostly by 
the 2-10 keV X-rays. Therefore, an increase in fn for a given M achieves: (i) an increase in 
the plasma column, density and luminosity and (ii) a relative decrease of the ionizing hard 
X-ray flux (smaller ct ox ). With these changes over FKCB, in the present treatment iron is 
not fully ionized even at the smallest radii (r ~ 10r s ), leading to w(Fexxv) ~ 0.7c. This 
partial ionization of the plasma, coupled with the increased column, reduces the ionizing flux 
that reaches further out into the wind, especially for high 9, leading to a bootstrap of less 
ionization and increasingly higher soft X-ray opacity. Then, because of the ensuing severe 
reduction of ionizing photons in the E ~ 0.1 — 2 keV range, the Civ ions form at smaller 
radii (higher velocities), so that the r~ 2 increase of the photon flux with E > 64 eV (the 
ionization potential of C iv) off-sets the photon depletion due to photoelectric absorption by 
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the partially ionized plasma. 

Recent spectroscopic studies of BAL QSOs have indicated likely correlations between the 
maximum outflow velocity of the UV/X-ray absorbers and the spectral indices (G06;C09). 
For comparison we show in Figure |3^l the modeled outflow velocities for C IV and Fe xxv with 
different values of (r, a ox ) for rh = 0.5 and 9 = 50°. It is seen for a ox = —2 that Fexxv ve- 
locities correlate strongly with F (solid curves) allowing for velocities in the range 0.3 < 
y(Fexxv)/c < 0.8 consistent with the X-ray outflow velocity observations in APM 08279+5255 
(C09), while the Civ velocity is virtually unaffected. This is because for the steeper X-ray 
spectra fewer ionizing photons are available to produce highly-ionized species (e.g. Fexxv of 
ionization potential ~ 9 keV) and the relevant ions are found at smaller distances (higher 
velocities) than in the case of harder spectra. This does not affect significantly the overall 
ionization of the wind leaving the C IV transition at roughly the same distance. However, a 
change in a ox , affects strongly the maximum column position of both Fe xxv and C IV, as de- 
scribed above. For constant T(= 2), the velocities of both these transitions correlate strongly 
with a ox (dashed curves), ranging between 0.1 < t>(Fexxv)/c < 0.8 and 0.0! < v(Cw)/c 



<^ 0. 15 for —2.1 <^ a n x <> —1.6 qualitatively consistent with UV data (ILaor fe Brandtl 12002 



G06; iFan et al.l 120091 ) . Radiation forces are often invoked to explain these correlations (with 
X-ray shielding necessary for high velocities, e.g. MCGV;PSK). In contrast to these models, 
high-velocity flows are always present in our model, but only the steep a ox (<^ —2) allows the 
relevant ions (e.g., Fexxv and Civ) to form in their small-r, high-v regions which are oth- 
erwise overionized (c.f. FKCB). While some (narr ow- line) Seyferts with — 1.6 < a m ^ — 1. 1 



appear to exhibit X-ray outflows with v/c < 0.15 (IDadina et al.ll2005t iTombesi et al.1120 10'). 



these are systematically slower [an d generally substantially slower (IHolczer. Behar fe Kaspi 



20071 : IHolczer. Behar fc Aravil2010h ] than those of the BAL QSOs with a ox £ -1.6. Thus, 
we propose a ox as the defining parameter that determines the velocities of the UV/X-ray 
absorption features in AGNs. 

While the intrinsic MHD wind ionization structure is determined by a ox , the observables, 
i.e. the velocity widths/shifts of the ions depend strongly also on the observers' inclination 
angle 9. In Figure [3b we present the LoS velocity of the Civ and Fexxv transitions for 
various 9 with V = —a ox = 2. Because of the specific geometric shape of the magnetic field 
lines and ionization equilibria, characteristic ion velocities vary for different LoS angles [see 
also Fig. la for their positional transitions along various LoS from 9 = 80° (innermost) to 
30° (outermost)]. In this fiducial model we find that w max (Fexxv) ~ 0.6c at 9 ~ 50° and 
^max(Clv) ~ 0.15c at 9 ~ 70°. At larger angles the velocities decrease but the integrated 
columns are so high that it is doubtful these features are observable. Similar diagrams can 
be computed for different values of the parameters (rh, T, a ox ) and can be directly compared 
to observations to assess the fundamental assumptions of these models. 
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4. Summary & Discussion 



We have demonstrated that purely MHD disk- winds with n oc r" 1 , originally proposed to 
account for the X-ray AMDs in Seyferts, can also encompass combined high-velocity UV/X- 
ray absorber properties as diverse as those of BAL QSOs, with those of APM 08279+5255 as 
a template. This is extremely important in view of the winds' scale invariance, with the qual- 
itative differences in the absorber properties between Seyferts and QSOs attributed mainly 
to their different m and a m . Given the well documented correlation of AGN UV- lumi nosity 
(a proxy for rh) with a ox ( jLaor fc Brandtl l2002t ISteffen et al.l 120051 ; Fan et al.l 120091 ). this 
model implies AGN absorber structure that depends essentially on a single parameter a ox . 
However, the observables, e.g., columns and velocities, depend additionally on the LoS angle 
9, reproducing the QSO BALs Civ (UV) and Fexxv (X-ray) properties only for sufficiently 
large 9, as usually considered. 

Our calculations show that Civ forms at r/r s ~ 200 — 700 with corresponding velocity 
v (Civ) <, 0.1c. This value, along with the Fexxv velocity v (Fexxv) ~ 0.6c at r/r s ~ 5-40, 
are consistent with those observed in APM 08279+5255 (C09). The shielding of the plasma 
from the X-rays at r > 100r s , necessary to produce the high velocity Civ absorption in 
radiation-driven wind models (MCGV;PSK), is in our case naturally provided by the faster 
components of the same wind launched from even smaller radii. Ionization equilibrium as a 
result of steeper a ox allows the relevant charge states to form closer to the central engines 
where the wind is faster. 

For simplicity, we have ignored he re a number of physical processes, e.g. radiation pres- 
sure (see PSK) and thermal instability ( iKrolik. McKee fc Tarterlll98ll : iHolczer. Behar fc Kaspi 
20071). likely to h ave an impact on our MHD-wind properties, that need to be implemented 
(e.g. |Progall2003l ). However, the broader validity of our models, gauged by the AMD depen- 
dence on £, will be decided by observations and quantitative analysis such as those presented 
in B09. It is also encouraging that the ionization properties of certain X-ray absorbers are 
consistent with this picture an d, in fact, magnet ic-driving of disk- wi nds has been favorably 
argue d for GRO J1655-40 (e.g. lMiller et al.ll2008h and NGC 4151 (e.g. ICrenshaw k Kraemer 
20071 ). for example. We anticipate the upcoming A stro-H mission to contribute significantly 
to this goal by providing more detail on the Fe-K component of the wind, and thus to further 
clarify our picture of AGN structure. 
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- 10 - 



REFERENCES 

Behar, E., 2009, ApJ, 703, 1346 (B09) 
Blandford, R. D & Payne, D. G. 1982, MNRAS, 199, 883 
Brandt, W. N., Laor, A., & Wills, B. J. 2000, ApJ, 528, 637 
Castor, J. I., Abbott, D. C. & Klein, R. I. 1975, ApJ, 195, 157 

Chartas, G., Brandt, W. N., Gallagher, S. C, & Garmire, G. P. 2002, ApJ, 579, 169 

Chartas, G., Brandt, W. N., Gallagher, S. C. 2003, ApJ, 595, 85 

Chartas, C, Eracleous, M., Dai, X., Agol, E., & Gallagher, S. C. 2007, ApJ, 661, 678 

Chartas, G., Saez, C, Brandt, W. N., Giustini, M., & Garmire, G. P. 2009, ApJ, 706, 644 
(C09) 

Contopoulos, J., & Lovelace, R. V. E. 1994, ApJ, 429, 139 (CL94) 
Crenshaw, D. M. et al. 2003, Annu. Rev. Astro. Astrophyps., 41, 117 
Crenshaw, D. M. & Kraemer, S. B. 2007, ApJ, 659, 250 

Dadina, M., Cappi, M., Malaguti, G., Ponti, C, & de Rosa, A. 2005, A&A, 442, 461 

Dai, X., Kochanek, C. S., Chartas, G., Kozlowski, S., Morgan, C. W., Garmire, G., & Agol, 
E. 2010, ApJ, 709, 278 

Everett, J. E. 2005, ApJ, 631, 689 

Fan, L. L., Wang, H. Y., Wang, T., Wang, J., Dong, X., Zhang, K., Cheng, F. 2009, ApJ, 
690, 1006 

Fukumura, K., Kazanas, D., Contopoulos, I., & Behar, E. 2010, ApJ, 715, 636 (FKCB) 

Gallagher, S. C, Brandt, W. N., Sambruna, R. M., Mathur, S., & Yamasaki, N. 1999, ApJ, 
519, 544 

Gallagher, S. C, et al. 2006, ApJ, 644, 709 (G06) 
Green, P. J., & Mathur, S. 1996, ApJ, 462, 637 
Holczer, T., Behar, E. & Arav, N. 2010, ApJ, 708, 981 
Holczer, T., Behar, E., & Kaspi, S. 2007, ApJ, 663, 799 



- 11 - 



Hewett, P. C., & Foltz, C. B. 2003, A J, 125, 1784 
Kallman, T., & Bautista, M. 2001, ApJS, 133, 221 
Krolik, J., McKee, C, & Tarter, 1981, ApJ, 249, 422 

Kopko, M., Turnshek, D. A., & Espey, B. R. 1994, in IAU Symp. 159, Multiwavelength 
Continuum Emission of AGN, ed. T. Courvoisier & A. Blecha (Dordrecht:Kluwer), 
450 

Laor, A. & Brandt, W. N. 2002, ApJ, 569, L641 

Miller, J. M., Raymond, J., Reynolds, C. S., Fabian, A. C, Kallman, T. R., & Homan, J. 
2008, ApJ, 680, 1359 

Murray, N., Chiang, J., Grossman, S. A. & Voit, G. M. 1995, ApJ, 451, 498 (MCGV) 

Pounds, K. A. et al. 2003, MNRAS, 345, 705 

Proga, D., Stone, J. M., & Kallman, T. R. 2000, ApJ, 543, 686 

Proga, D. 2003, ApJ, 585, 406 

Reeves, J. N., O'Brien, P. T., & Ward, M. J. 2003, ApJ, 593, L65 
Reeves, J. N. et al. 2009, ApJ, 701, 493 
Sim, S. A. 2005, MNRAS, 356, 531 

Sim, S. A., Long, K. S., Miller, L., & Turner, T. J. 2008, MNRAS, 388, 611 
Steffen, A. T. et al. 2005, ApJ, 131, 2826 

Sim, S. A., Miller, L., Long, K. S., Turner, T. J., & Reeves, J. N. 2010, 404, 1369 

Srianand, R. & Petitjean, P. 2000, A&A, 357, 414 

Tananbaum, FL, et al. 1979, ApJ, 234, L9 

Tombesi, F. et al. 2010, accepted to A&A, arXiv: 1006.28581 

Weymann, R. J., Morris, S. L., Foltz, C. B., & Hewett, P. C. 1991, ApJ, 373, 23 



This preprint was prepared with the A AS IATgX macros v5.2. 



